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Abstract 

The problem of spin precession in a time-dependent magnetic field is con- 
sidered in the adiabatic approximation where the field direction or the an- 
gular velocity of its rotation is changing slowly. The precession angles are 
given by integrals in a way similar to the semi-classical approximation for the 
Schrodinger equation. 
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Resume. Le probleme de la precession du spin dans un champ magnetique 
variable est envisage selon 1' approximation adiabatique dans laquelle la direc- 
tion du champ ou la vitesse de sa rotation varie lentement. Les angles de la 
precession sont donnes par quadratures de maniere similaire a l'approximation 
semi-classique de l'equation de Schrodinger. 
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Version francaise abregee 

Le probleme de la precession du spin dans un champ magnetique variable a ete 
resoiu explicitement pour deux cas particuliers importants (Landau and Lifshitz 
1974, Sect. 114 ): (i) le champ magnetique ne change pas de direction, (ii) le 
champ de valeur constante tourne autour d'un axe donne a une vitesse angulaire 
constante. En outre, des solution analytiques sont obtenues pour certaines vari- 
ations du champ (cf. Rosenfeld et al. 1996). La presente etude a pour but de 
presenter une solution approchee de l'equation de Bloch (1) decrivant la precession 
du spin. L'approximation est valabie si la direction du champ ou de sa vitesse change 
adiabatiquement; le resultat principal est enonce par l'equation (16). 
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Nous considerons l'equation de Bloch dans le cas ou la dependance du champ 
magnetique B(t) de temps t est arbitraire. La solution fondamentale est donnee 
par une matrice unitaire U(t,to) satisfaisant l'equation (3). Cette matrice peut 
etre ecrite en termes de deux fonctions complexes (£,r)), qui satisfont aux equations 
differentielles (5). Des equations de ce type et des solutions approchees ont ete 
considerees precedemment. Le probleme peut etre reduit a celui d'un oscillateur a 
frequence variable et complexe, eq. (7). II est a noter que le carre de la frequence 
d'oscillation dans l'eq. (7) est une constante reele positive pour les cas simples, mais 
s'avere etre generalement complexe. C'est pourquoi la methode traditionelle WKBJ 
ne peut etre appliquee directement. 

Dans le but de trouver une bonne approximation, la matrice U dans l'eq. (9) 
est decrite par une rotation d'un angle 2a autour de la vecteur unitaire n = 
(sin 7 cos 5, sin 7 sin 5, cos 7). Le probleme consiste en recherche de la dependance 
temporelle des trois angles (a, 7, 5) a partir des equations differentielles (11)-(13). 

L'hypothese cruciale, qui simplifie le probleme, est que revolution de 7 est tres 
lente. En d'autres termes, nous supposons que le vecteur n varie adiabatiquement. 
Les equations resultantes sont integrees immediatement par quadratures, eq. (16). 
II est remarquable que pour les deux cas mentionnes dans le premier paragraphe, 
1' approximation est exacte, puisque soit (i) Bj_ = 0, l'axe de rotation est toujours 
aligne sur le champ, soit (ii) (p, B z et Bj_ sont fixes, de sorte que Tangle entre le 
champ et l'axe est preserve, ainsi que Tangle 7. 

L' approximation edt effective sous deux conditions alternatives: soit (i) la di- 
rection du champ B, ou (ii) la direction du vecteur (B x B), varient lentement 
autour d'un vecteur e. Dans les deux cas, le vecteur e indique l'axe z pour le repere 
approprie et pour le choix des parametres de rotation dans l'equation (9). 

La forme de la solution (16) obtenue ressemble a T approximation semi-classique 
(ou eikonale) appliquee a la diffusion a haute energie, pour laquelle le potentiel de 
diffusion etait suppose etre une fonction d'espace a variation lente. 

* * 
* 

The problem of spin precession in variable magnetic field is fundamental for the 
theory of nuclear magnetic resonance (Abragam 1961). Its explicit solution is known 
for two important particular cases (Landau and Lifshitz 1974, Sect. 114 ): i) the 
magnetic field does not change its direction, ii) the field has a constant magnitude 
and is precessing around a given axis with a constant angular velocity. Besides, 
analytical solutions have been obtained for a number of special field shapes (Rosen 
and Zener 1932, Bambini and Berman 1981, Silver et al. 1985, Rosenfeld et al. 
1996). 

As soon as one deals with the Cauchy problem for two coupled linear first-order 
differential equations, or equivalently, with a second-order equation, an adiabatic 
approximation may be employed, similarly to the standard semi-classical approach 
to the Schrodinger equation. An adiabatic approximation for the Bloch equation 
was the purpose of the present work, and the result is presented in Eq.(16). The 
key to the desired approximation is a proper group-theoretical parametrization in 
Eqs.(9-10), which enables one to handle complexities specific to the problem. 
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The spinor wave function of a neutral spin-^ particle satisfies the Bloch equation, 

id^/dt = (B • a)i/j, (1) 

where B = (B x , B y , B z ) = fi>B(t), B(t) is the variable magnetic field vector, [i is the 
particle magnetic moment, and a are the Pauli matrices. The fundamental solution 
of Eq. (1) is given by a unitary 2x2 matrix U, 

m = u(t,toMt ), (2) 

idU/dt = (B • a)U, U(t , t ) = /• (3) 

(Here / is the unit matrix.) Two complex functions (£,77) are used to represent U, 

£ e -H3 ^-03 \ _ r t 

I to 



In view of Eq. (3), these functions satisfy the following equations, 

i£ = -bn, iv = -k, |e| 2 + H 2 = i- (5) 

where 

b = (B x - iB y )e 2lP = B ± e 2il3 - 2i(p , B\ = B 2 - B 2 Z , tan 2 V = B y /B x . (6) 

Equations of this type and various approximations were considered previously (e.g. 
by Popov 1962, Marinov and Popov 1979, Rosenfeld et al. 1996). The problem may 
be further reduced to that for an oscillator with a variable (complex) frequency, 

f + n 2 (t)f = o, f = b-v% (7) 

where 

n 2 = B]_ + uj 2 + id), 
u; = -tb/2b = B z -<p- — log Si, (8) 

i.e. 2ip = (B x B y — B x B y )i ' B\ is the angular velocity of the field rotation around 
the z-axis. Note that Q 2 is a real positive constant for the cases presented in the 
textbook, but is complex in general. Therefore the standard WKBJ method cannot 
be applied to Eq. (7) directly. 

A consistent adiabatic approximation for the problem concerned was proposed 
by Bender and Papanicolaou 1988 and further developed by Papanicolaou 1988. It 
was assumed that the field time dependence was slow, i.e. the specific variation time 
T was large, and an expansion in powers of 1/T was constructed. An alternative 
approach, based upon the Magnus expansion, was presented by Klarsfeld and Oteo 
1992. Unlike the above mentioned (and a number of other) works, we do not assume 
a slow time dependence of the field. In particular, the result given below is exact for 
the field rotation with an arbitrarily high (but stable) angular velocity. Actually, we 
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employ geometrical arguments, which were previously developed for the adiabatic 
approximation by Nakagawa 1987. With all that, to the best of our knowledge, the 
result given in Eqs. (9), (16) was never published before. 

The complex equations (5) are equivalent to real equations for angles which 
appear when U is treated as the spinor representation of the rotation group. Let us 
write U as a product of two matrices representing rotations, 

U = R z (ai)R n (a), R n = Icosa — i(n ■ a) sina. (9) 

Here matrix R z {a\) has the diagonal elements (e iai , e~ iai ) only and represents the 
rotation around the z-axis, and R n (a) represents a rotation by an angle 2a around 
the direction of a unit vector n = (sin 7 cos 5, sin 7 sin 5, cos 7) . The factor containing 
a±(t) transforms the system to a rotating frame. The time dependence of a± is chosen 
in such a way as to make the variation of n slow. Setting 

ai = -<p + 6/2, (10) 

we get the following differential equations for the angles (a, 7, 5), 

d sin7 = z/sin(27 — 70), (11) 

7 tan a = 2z/sin(7 — 7), (12) 

<5sin7 = 2t/sin(7o — 7), (13) 

where the new notations are 



z/(t) = v / V 9 2 -2 v 3 J B 2 + B 2 , cot>y (t) = (<p-B x )/B ± . (14) 
The initial conditions corresponding to (3) are 

a(t o ) = 0, 7(fo)=7o(*o), S(t ) = <p(t ). (15) 

The starting precession velocity of n, as given by Eqs. (12), (13), is vanishing, and 
the crucial assumption is that it is small for all t. In the other words, we assume 
that the vector n is moving adiabatically. In contrast to the standard adiabatic 
approximation, n does not coincide with the field direction if its change is rapid. As 
seen from Eq. (14), n takes a position intermediate between B and its instantaneous 
rotation axis (B x B). The approximate solution of (11)-(13) is 

a(t, t ) = f v(r)dT, 7 (t, t ) = 7o (f), S(t, t ) = <p(t ). (16) 

Jto 

For two cases mentioned in the first paragraph the approximation is exact, since 
either i) B± = 0, so 70 = 0, or ii) ip, B z , and B± are fixed, so 70 is a constant too. 

In order to estimate the accuracy of the approximation, let us assume that 70 is 
small and 7 ~ 70 . Respectively, equation (12) is linearized and solved explicitly 



7 PS 2(7 — 7)0; cot a, 



7(Mo) ~ 7o(*) + 



to 



7o(r) sin 2 a(r, t )dr 



/sina 2 (t,t ), (17) 
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where a is given by (16). The approximation is effective under two alternative 
conditions: either (i) the field direction, or (ii) the direction of the vector (B x B), 
is changed slowly around a fixed vector e. In both cases, the vector e indicates the z 
axis for the appropriate coordinate system and the choice of the rotation parameters 
in Eq. (9). A remarkable case is that of a fastly rotating field, where (p ^> B±, and 
| cot 7o| 3> 1. As a result, the deviation of n from the z direction is small, even if the 
field rotation is not stationary. The obtained form of the solution (16) resembles the 
semi-classical and the eikonal approximation to the high-energy scattering (Landau 
and Lifshitz, 1974, Sect. 131), since it is valid if z> u 2 . 

The validity of our approximation was checked by comparing it with the exact 
result for a field where the analytical solution was obtained by Rosen and Zener 
1932. That field has a constant x-component and a pulse in the ^-direction, 

B «>=H 8 °-;^W o )- < i8 > 

where (3q and ( are constant dimensionless parameters. Assuming that the initial 
state was polarized to the x-axis, the spin-flip probability W was calculated in the 
time asymptotics, i.e. t, \to\ ^> T. In the limit of (3q — > oo, corresponding to a large 
T for a fixed background field, the exact and the approximate solutions coincide, 

W = tVf3l + 0{{3^). (19) 

Moreover, the numerical calculations showed that the relative error of our approxi- 
mation in W is less than 1% for all values of @ > 1-5 and ( > 1. 
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